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Phase diagram of a bosonic ladder with two coupled chains. 



OO 
O 

o 

(N 

X> 

<D 

IX. 



O 

B 

o3 



c 

o 
o 



> 

On 

o 

O 

00 

O 



X 



Meetu Sethi LuthraB and Tapan MishrE0 
Indian Institute of Astrophysics, II Block, Kormangala, Bangalore, 560 034, India. 

Ramesh V. Pa§ 

Department of Physics, Goa University, Taleigao Plateau, Goa 403 206, India. 

B. P. Dasl 

Indian Institute of Astrophysics, II Block, Kormangala, Bangalore, 560 034, India. 

(Dated: February 5, 2008) 

We study a bosonic ladder with two coupled chains using the finite size density matrix renor- 
malisation group method. We show that in a commensurate bosonic ladder the critical on-site 
interaction (Uc) for the superfluid to Mott insulator transition becomes larger as the inter-chain 
hopping (f^)increases. We analyze this quantum phase transition and obtain the phase diagram in 
the t± — U plane. 
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I. INTRODUCTION 



Quantum phase transitions in ultracold atoms provide 
important insights into the behaviour of matter at very 
low temperatures [H, 0, S 0] ■ An important example of 
this class of transitions is the transition from a super- 
fluid(SF) to a Mott insulator (MI) which has been ob- 
served in cold bosonic atoms in 3D optical lattices [f| Q 
as predicted by Jaksch et al [7[. Subsequently, this tran- 
sition has been observed in a ID optical lattice [||. De- 
tailed theoretical studies of this transition have been car- 
ried out by Pai et al [a l2(- An important question to 
address is how do the characteristics of this transition 
alter in going from one to two dimensions. In coupled 
bosonic chains, competition exists between the ratio of 
atomic interactions to the intra-chain hopping and the 
inter-chain hopping. A large value of the former favors 
a Mott insulator state overcoming the effect of the inter- 
chain hopping while if the latter dominates it would tend 
to delocalize the bosons and drive the system to a super- 
fluid state. It is not practical to perform numerical stud- 
ies for the above mentioned transitions in a very large 
number of coupled chains, so one must restrict to a fi- 
nite number of such chains. The aim of the present work 
is to study the effect of the inter-chain hopping on the 
SF-MI transition for a bosonic ladder consisting of two 
coupled chains. Although a substantial amount of theo- 
retical and numerical work has been done in this direction 
for the case of spin-less fermionic ladders and spin lad- 
ders [1(1 HH, EH ,no work has been done to our knowledge 
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for the bosonic ladders except a recent work using the 
Bosonization method. [HI, The Hamiltonian of the 
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FIG. 1: Schematic picture of a two-leg bosonic ladder, t and 
t± are, respectively, inter-chain and intra-chain hopping am- 
plitudes. 

bosonic ladder(as shown in Fig[T])is given by 
H = -*y^(aj, a ffl <+i,ot + ft- c ) 

u 



i,a> 

-t± ^(a!,i fl i,2 + h.c). 



(i) 



In this model (HJ a\ a (aj, a ) represents bosonic creation 
(annihilation) operator for the site i of the chain with 
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index a = 1, 2. t and {/ are the intra-chain hopping am- 
plitude between the nearest neighboring sites of chain a 
and the on-site interaction between the bosons respec- 
tively. The last term in this model |T]) represents inter- 
chain hopping with an amplitude t±_ between correspond- 
ing sites on the two chains. We set our energy scale by 
taking t = 1. 

This model has been studied using the Bosonization 
technique [H, [3] at or close to commensurate filling of 
one boson per site. This study predicts a transition from 
a Mott insulator to a superfluid phase when the inter- 
chain hopping is increased and it is in the Beresenskii- 
Kosterlitz-Thouless (BKT) universality class at commen- 
surate filling. In the present work we verify these predic- 
tions and thus complement the earlier analytical results. 
For this purpose, we study the variation of critical on-site 
interaction Uc for the superfluid phase to the Mott insu- 
lator phase transition with the change in the inter-chain 
hopping amplitude ij_ using the finite size density matrix 
renormalization group method (FSDMRG) [IEEE!] and 
obtain the phase diagram in the (t± — U) plane. To the 
best of our knowledge, the present work is the first ap- 
plication of FSDMRG method to bosonic ladders. 

The remaining part of the paper is organized in the 
following manner. The FSDMRG method in the context 
of bosonic ladders is briefly described in section [TXJ Our 
results are described and discussed in Section HTT1 and our 
conclusions are stated in Section ITVl 



sufficient for the values of U considered here The 
value of M is chosen such that the truncation error in 
our calculation is always less than KG 5 . 
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FIG. 2: A scheme of superblock configuration for the FS- 
DMRG algorithm 



II. FSDMRG METHOD 

We use the FSDMRG technique to obtain the ener- 
gies and the correlation functions of the ground state. 
This method is very efficient and has proven to give accu- 
rate results for ID quantum lattice systems and has been 
applied to low-dimensional strongly correlated fermionic 
and bosonic systems We give below some perti- 

nent details of this method adapted to the system of two 
coupled chains that we have considered. 

We begin with a super-block configuration 
B l L • • B r L of L rungs as shown in FigJ2] The 

2 2 

left B l L and the right block B r L have 4 — 1 rungs 

each and the • represents one rung of two sites, one 
from each chain. Thus in every iteration, the new left 
and right blocks are B 1 ^ = B\ • and B 1 ^ = • B r I± _ 

"2 ~2~ 1 _ "2" 

respectively. We increase the size of the system in 
every iteration by adding two rungs which increases 
the number of lattice sites by 4. To keep the density 
p = 1 fixed, we also increase the number of bosons in 
the system by 4. The truncation of states of left (right) 
block in each iterations corresponds to choosing AI 
highest weighted states out of 2 x nmax X M of the 
left (right) density matrix. Here nmax is the number of 
states kept at each site, which is in general infinity, but 
we truncate it for a feasible numerical calculation. Wc 
keep n m ax = 4 in this calculation which is found to be 



III. RESULTS AND DISCUSSION 

The Mott insulator phase has a finite gap in its energy 
spectrum. The single particle gap is defined as 

G L = E L (N+1)-E L (N)-(E L (N)-E L (N-1)) (2) 

where E^{N) is the ground-state energy of two leg 
bosonic ladder with length L having N bosons. The MI 
phase is signaled by the opening up of the gap G^^oo. 
However, Gl is finite for finite systems and we must ex- 
trapolate to the L — > oo limit, which is best done by 
using finite-size scaling Q. In the critical region, i.e., SF 
region, the gap 

GzwL-V^/O. (3) 

where the scaling function f(x) ~ x, x — > and £ is the 
correlation length. £ — > oo in the SF region. Thus plots 
of LGl versus U, for different system sizes L, consist of 
curves that intersect at the critical point at which the 
correlation length for L = oo diverges and gap Goo van- 
ishes. The phase diagram, as discussed below, is obtained 
from such plots. 

It is now well known that the single chain Bose- 
Hubbard model with density p = 1 shows a SF-MI tran- 
sition with the critical on-site interaction Uc ~ 3.40. 
In order to understand the effect of the inter-chain hop- 
ping on this transition, we varied t± from to 20 and 
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obtained the corresponding critical on-site interaction 
U c (t±) for the SF-MI transition. We found that U c 
increases with t± and saturates in the limit t±_ — > oo. 
These results are highlighted by the plots of scaling of 
gap LGl versus U for different values of t± and lengths 
L. For example, in the Fig|3] we plot LGl versus U 
for t± = 0.4. The coalescence of LGl curves for dif- 
ferent values of L below U < 6.6 demonstrates the SF- 
MI transition with Uc(t± = 0.4) ~ 6.6 which is much 
larger than the corresponding value for the single chain 
Uc(t± — 0) ~ 3.4. FigHI represents similar plots for 
t± = 1. For this case the critical on-site interaction in- 
creases further to U~c{{t± = 1) ~ 7.9. 
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FIG. 3: Scaling of gap LGl as a function of U for t± = 0.4 
and different lengths. The coalescence of curves for different 
lengths for U < Uc ~ 6.6 shows a superfluid phase and a 
Mott insulator with finite gap for U > Uc- 

From similar plots of LGl versus U, we obtain the 
phase diagram for model |(TJ) in the t± — U plane and 
it can be seen in FigJSJ Uc for the SF-MI transition 
initially increases sharply as the inter-chain hopping t± 
increases. This phase diagram verifies the prediction of 
MI-SF transition with respect to increase in the inter- 
chain hopping t± [l(| • For higher values of t± , Uc tends 
to saturate. For t± ^> t, U, each rung has two one particle 
states: corresponding to bonding or anti-bonding. As 
predicted in the Bosonization [10| study, this problem 
then maps onto a single chain Bose Hubbard model with 
commensurate density p = 2 and on-site interaction U jl. 
To confirm this prediction we plot the variation of Uc 
with respect to t\_ in Fig[Sl Critical on-site interaction 
Uc for large t± converges to a value equal to 12.5 ± 0.3. 
Plotting LGl versus U for single chain Bose-Hubbard 
model for p = 2 in FigJ7]we find that Uc ~ 6.3, which 
is one half the converged value of Uc(t± = oo) ~ 12.5 
for the bosonic ladder confirming the prediction made in 
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FIG. 4: Scaling of gap LGl as a function of U for t± — t — 1.0 
and different lengths. The coalescence of curves for different 
lengths for U < Uc ~ 7.9 shows a superfluid phase and a 
Mott insulator with finite gap for U > Uc- Comparing this 
figure with Fig(3] we observe that the critical Uc increases 
with tj_. 

Bosonization study [Iol]. 

The Mott insulator to superfluid transition is found 
to be Beresinskii-Kosterlitz-Thouless (KT) universality 
class at commensurate filling The correlation func- 
tion that characterizes the superfluid phase is given by 
r a (r) = {a[ a a,i+r,a) which decays as a power law in the 
limit t — > oo. Here the expectation value is taken with re- 
spect to the ground state. However, in the Mott insulator 
phase it has an exponential decay due to the finite gap in 
the energy spectrum. The power law decay of this corre- 
lation function has been obtained using the Bosonization 
method [l(| and it is predicted to go as 

with the Luttinger Liquid parameter K s stated to be 1 
at the superfluid to Mott insulator transition point. 

In order to obtain the Luttinger Liquid parameter K s 
we fix U = 6 and vary the inter-chain hopping tj_ and 
obtain the Mott insulator to superfluid transition. The 
scaling of gap LGl as a function of t± is given in FigEl 
The critical inter-chain hopping t± = 0.24 ± 0.05 for the 
MI to SF transition. The correlation functions T a (r) for 
C/ = 6 different values of t± are given in FigEJ The Lut- 
tinger Liquid parameter K s which is obtained by fitting 
T a (r) with the expression given in EqH] is plotted as a 
function of i_L in FigQIj] From these values the critical tj[_ 
for which K s = 1 is given by 0.3 ±0.03 which is consistent 
with values obtained from the scaling of the gap. 

The results we have obtained could have experimen- 
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FIG. 5: Phase diagram of model ([T| as a function of inter- 
chain hopping t± and on-site interaction U for density p = 1. 
Note that we have set intra-chain hopping t = 1. 



FIG. 7: Scaling of gap LGl as a function of U for the single 
chain Bose-Hubbard model with density p = 2. The coales- 
cence of curves for different lengths for U < Uc ~ 6.3 shows 
a superfluid to a Mott insulator transition. 




FIG. 6: Variation of critical on-site interaction Uc with re- 
spect to inter-chain hopping t±. Uc increases sharply for 
small values of t± and saturate to 12.5 ± 0.3 as t± — » oo. 



FIG. 8: Scaling of LGl as a function of t±_ for U = 6. 
IV. CONCLUSIONS 



tal implications. It is now possible to prepare bosonic 
ladders by growing optical superlattices in the form of 
double well potential along one direct ion [l8j. The tun- 
neling between the double well potential will control the 
inter-chain hopping. By changing dynamically the opti- 
cal lattice parameters one can control all the interaction 
and hopping parameters of the model (|TJ) . 



We have studied the ground state properties of a two 
chain bosonic ladder with commensurate filling of one 
boson per site using the finite size density matrix renor- 
malization group method. The critical on-site interaction 
for SF-MI phase transition increases sharply for small 
values of inter-chain hoping amplitude t±. However, it 
saturates in the limit t± — > oo. Thus in the presence of 
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FIG. 9: Power law decay of T a (r) for various values of t± for 
U = 6. 
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FIG. 10: Variation of Luttinger liquid parameter K 3 as a 
function of t± near MI to SF transition for (7 = 6. 
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